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We present an attempt of global analysis of Semi-Inclusive Deep Inelastic Scattering (SIDIS) 
lp^ — > £'irX data on single spin asymmetries and data on left-right asymmetry An in p^p —¥ ttX in 
order to simultaneously extract information on Sivers function and twist-three quark-gluon Efremov- 
Teryaev-Qiu-Sterman (ETQS) function. We explore different possibilities such as node of Sivers 
function in x or k± in order to explain "sign mismatch" between these functions. We show that 7T 
SIDIS data and tt° STAR data can be well described in a combined TMD and twist-3 fit, however 
" n ± BRAHMS data are not described in a satisfactory way. This leaves open a question to the 

solution of the "sign mismatch" . Possible explanations are then discussed. 
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03 . I. INTRODUCTION 
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Single transverse spin asymmetries (SSAs) are a rich source of information on internal partonic structure of the 
nucleon The exploration of the underlying mechanisms has led us to realize that the SSAs are sensitive probes 
of the parton's transverse motion. There are two different yet related QCD factorization formalisms to incorparate 
i-C 1 such transverse components of the parton's momentum and to describe the observed asymmetries: the transverse 
^4 momentum dependent (TMD) factorization and the collinear twist-three factorization approaches. 

For processes such as single inclusive hadron production in proton-proton collisions, p' p — > hX , which exhibits only 
one characteristic hard scale, the transverse momentum P 2 ± ^> Aq CD of the produced hadron, one could describe the 
SSAs in terms of twist-three quark-gluon correlation functions 0t3| • One of the well-known examples is the so-called 
Efremov-Teryaev-Qiu-Sterman (ETQS) function. Phenomenological extractions were performed in different papers 
[1,0. On the other hand, for processes such as Semi- Inclusive Deep Inelastic Scattering (SIDIS) which possesses 
^ | two characteristic scales, photon's virtuality Q and Ph± of the produced hadron, one can use a TMD factorization 
formalism [H ED] 

in the region Aq CD < P 2 ± <C Q 2 and describe asymmetries with TMD functions. One of the most 

important TMDs is the Sivers function f^f [l2|, [l3| which describes sin^^ — </> s ) modulation in SIDIS on transversely 
l/"} ' polarized target [3]. Sivers functions have been extracted from SIDIS experimental data by various groups [IB - fl9l |. 

These two formalisms are closely related to each other, and have been shown to be equivalent in the overlap region 
where both can apply poU^ . The relevant functions - the Sivers function and the ETQS twist-three function - are 
| connected through the following relation (23l-|25jp1: 



X 



T q , F (x,x) =-[ d 2 k ± ^f^(x,kl)\ smis , (1) 

where the subscript "SIDIS" emphasizes that the Sivers function is probed in the SIDIS process. Color gauge invariant 
nature of TMDs manifests itself in the fact that TMDs are process dependent and an important consequence of this 
process-dependence is a prediction [27[ : 

/itIsidis = —firloY ■ (2) 

i.e., the Sivers function measured in SIDIS and Drell-Yan (DY) processes are exactly opposite to each other. Experi- 
ments are actively planning to measure and verify such a prediction. Somepreliminary phenomenological estimates 
on the SSAs of DY production [28 30] and solid theoretical developments [3l|, HH have been achieved. 

Recently it was found that left-hand-side (LHS) and right-hand-side (RHS) of Eq. (JTJ) have opposite signs if the 



corresponding functions are extracted from phenomenological studies of different experimental data [25[ , particularly 
the RHS f^j!(x, k 2 ^) from the SIDIS data, while LHS T q p(x, x) from pp data. We will refer this finding as "sign puzzle" 
or "sign mismatch". Whether it reflects the incompatibility of SIDIS and pp data within the current theoretical 
formalism, or reflects the inconsistency of our formalism itself, is a very important question and needs to be further 
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1 Note that this relation holds for the so-called unsubtracted part of TMDs, see Ref. |2(I . 
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explored both theoretically and experimentally. On the experimental side, the measurement for the SSAs of single 



inclusive jet and direct photon production [25[, and the single lepton production from a W^-boson decay [33|, |34| in pp 
collisions, the single inclusive jet and hadron production in Ip collisions (without identifying the final-state lepton) 
35 37} could be very helpful. The study of hadron distribution inside a jet could also be useful (38l - l4l| . 

In this paper we will try to make some first attempt on the theoretical (phenomenological) side: we will attempt 
to make a global fitting of both SIDIS and pp data with a more flexible functional form for the Sivers function (and 
the ETQS function), to see if we are able to describe all the data within our current theoretical formalism. Our naive 
starting point is based on the observation that the SIDIS and pp data typically covers slightly different kinematic 
region, in either momentum fraction x and/or the transverse components. Thus a sign-changing functional form in 
these kinematic space might be just needed to cure the "sign mismatch". One such possibility, e.g., a node in x 
region has already been indicated in Ref. [Hj]. We will consider SIDIS data from HERMES and COMPASS, and 
proton-proton data from STAR and BRAHMS. Let us emphasize that this is a first attempt to use simultaneously 
the TMD and collinear twist-three factorization formalisms in a global analysis of the spin asymmetry. 

The rest of our paper is organized as follows. In Sec. II we recall the basic formalisms needed to describe SIDIS 
data for semi-inclusive hadron production at low Ph j_ , and proton-proton data for inclusive hadron production at high 
Ph_L- In Sec. Ill we first introduce our more flexible parametrized functional form for the Sivers function, describe our 
fitting procedure. Particularly we explore the possibility of a node in x region, and investigate whether it can help 
resolve the "sign mismatch" problem. At the end of this section we briefly comment on the possibility of the node in 
kj_ space. We conclude our paper in Sec. IV. 

II. BASIC TMD AND COLLINEAR TWIST-3 FACTORIZATION FORMALISMS 

In this section we review the basic formulas for the spin asymmetries in both SIDIS and proton-proton processes. 
We start with the semi-inclusive hadron production at low Ph± in SIDIS, e(£) + A>(P, s±) —> e(£') + h(Ph) + X, which 
can be described by the TMD factorization formalism. The differential cross section for the so-called Sivers effect 
reads 0: 

cTo(x B ,y,Q 2 ) \Fuu,t + |s_l| sin(^ - s )F^" 0s) l , (3) 



dx B dy d(j> s dz h dcj) h Ph±dP h ± 

where cr = x ^ y Qi (l + (1 — y) 2 ) with a the fine structure constant, q — I — £' with q 2 = —Q 2 , and the usual SIDIS 
variables are defined as 

Q 2 P-q P-Ph (A , 

XB= 2p—q' V = —f Zft ~- () 

The Sivers asymmetry can be defined as the sin(0; l — <fi s ) modula: 

.sinffc-^,) „, . , u n Jd(phd<p s sm(<l)h-<j> s )(d(T((l>h,<i> s )-da((f>h,(f> s +n)) 
A U7 t v = 2(sm(</>, l - (f> s )) = 2+ r 12 12 u 72 — 7TTTTI — 2 — i — v\ > ( 5 ) 



where the subscript U stands for unpolarized lepton beam, and T for the transverse polarization of the target nucleon. 
In terms of structure functions one has 

a {XB,y,Q 2 ) Fuu.t 
The structure functions depend onig, Q 2 , %h an d P 2 ±> an d can be written as (43l |44| 

Fuu,t = x B ^2e 2 a / d 2 kj_ f a /A{xB,k 2 ± ) D h/a (z h , (P h x - z h kj_) 2 ) , (7) 
F^~^ =-x B J2 e lf d 2 k x *^f#(x B ,kl)D h/a (z h ,(P hX -z h k x ) 2 ) , (8) 

a 

where h = Ph±/\Phx\i fir ^ s the Sivers function, f a /A and D h / a are TMD parton distribution function (PDF) 
and fragmentation function (FF), respectively. All our definitions on the TMD functions and these expressions are 



consistent with the Trento convention 45j , which have been used in the experiments [46l . 147 
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On the other hand, for single inclusive hadron production at high Ph± in p^p collisions, A>(P, s±) + B(P') —> 
h(Ph) the spin-averaged cross section da = [da(s±) + da(— s±)]/2 is usually written in the collinear factorization 
formalism as, 



E h 



da ot s 
d^K = ~S 



f dz f dx' f dx 

/ —D h/C (z) / — f b/B ( x ') / —f a/A (x)H^ c (s,i,u)6(s + i + u), 



(9) 



where ^ a b c runs over all parton flavors, S = (P + P') 2 , f a /A( x ) and fb/B( x ') are the collinear PDFs, and D h / C (z) 
is the collinear FF. H a u b ^ fC are the well-known unpolarized hard-part functions for partonic scattering [HI, |4{|. s, i, 
and u are the usual partonic Mandelstam variables, for the final hadron of transverse momentum Pf l ± and rapidity y 
we obtain 



= xx'S, i — —xj ' zPh±\f~Se v , u — —x' / zPh±V~Se y . 



(10) 



The commonly used Feynman-x^ can be written as xf = smh(y). Note that the partonic x,x' and z are 

integrated over in Eq. (j9|. 

The spin-dependent cross section dAa(s±) = [da(s±) — da(— s ±)]/2 is given by the collinear twist-three factorization 
formalism: 



dAa(s ± ) 



d?P h 



Si vers 



^slPH^Yl j ^ D h/c( Z ) J ^rfb/B(x') J 



dx 

X 



1 

zu 



T a , F (x, x) - x—T a , F (x, x) 



H, 



Sivcrs 



(s, t, u)S (s + t + u) 



(11) 



where e a p is a two-dimensional anti-symmetric tensor with eia = 1 (and 621 = —1), T ay p(x,x) is twist-tree ETQS 
functions, and H^l™ r *(s 7 i, u) are the relevant hard-part functions which have been given in Refs. [||,[5(|- The subscript 
"Sivers" here is to remind that there are other types of contributions to the SSAs for the inclusive hadron production. 
What is written in Eq. (|TT|) is only the so-called soft gluon pole contribution [8|, and there could be soft fermion 
pole contribution [6[, and also the contribution from the twist-three fragmentation function [7| . Nevertheless, the 
extensive phenomenological study of the single inclusive hadron production has been performed for the soft gluon 
pole contribution [§, 9j, which indicates the soft fermion pole contribution is relatively small at least in the forward 
region where the asymmetry is the largest @ . Our study in the current paper will also concentrate on the soft gluon 
pole contribution, for which the relevant twist-three function - ETQS function T ay p(x,x) - has a close relation to the 
Sivers function as in Eq. (fTJ and thus we are able to perform a global analysis for both SIDIS and proton-proton data. 
We will comment on the contribution of the twist-three fragmentation function in the end of the next section. 
The SSA, Ajv, is given by the ratio of spin-dependent and spin-averaged cross sections 



A 



N 



E, 



dAa(s ± ) 
' d 3 P h 



da 
d?Ph' 



(12) 



The absolute sign of A n depends on the choice of frame and the coordinate system. In the center-of-mass frame of 
the incoming hadrons A and B, a convenient coordinate system (consistent with the experimental convention) is: the 
polarized nucleon A moves along +z, unpolarized B along — z, spin sj_ along y, and transverse momentum Pux along 
x-direction, respectively. In this frame 



(13) 



which should be used in Eq. pip . 



III. GLOBAL FIT OF THE SPIN ASYMMETRY: AN ATTEMPT 



So far all the phenomenological studies on the spin asymmetries in the market have been separated into two isolated 
parts. On one side, people use TMD factorization formalism to describe the SIDIS data for hadron production at 
low Phi., and they solely concentrate on SIDIS data, and do not include proton-proton data in the global fitting. 
The Sivers functions have been extracted as a result of such studies. On the other side, collinear twist-3 factorization 
formalism is used to describe the proton-proton data for single inclusive hadron production at high Ph±, and only 
proton-proton data are analyzed without inclusion of SIDIS data in the global fitting. The so-called ETQS functions 
have been extracted from such studies. However, as we emphasize in our introduction, Sivers and ETQS functions are 
closely related. Thus in this section we attempt to perform a global analysis of both SIDIS and proton-proton data 
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on the spin asymmetries. We will use TMD formalism to describe the SIDIS data in terms of the Sivers function. 
From the parametrization of the Sivers function, we obtain the functional form for ETQS function through Eq. ([1}. 
Then we use the collinear twist-3 formalism to describe proton-proton data in terms of our obtained ETQS function. 
In this way, we hope a single parameterization for the Sivers function could help us achieve a global fitting of both 
SIDIS and proton-proton data. We first introduce our parametrization for both Sivers and ETQS functions, then we 
present and discuss the results from our global fitting. We explore the possibility of node in x in details, and briefly 
comment on the possibility of node in fcj_ at the end of this section. 



A. Parametrization for the Sivers and ETQS function 

Following Refs. [13, EH, we parametrize both the spin-averaged PDF / m(x, fcj_) and FF D) l / a {z,p\,) with a 
Gaussian form for the transverse components: 



fa/A(x,k 2 ± ) = fa/A^^xe^/^, (14) 



D h/a (z,p 2 T ) = D h/a (z) -L-e-P?-/^), (15) 

such that they reduce to the usual collinear PDF f a /A( x ) an d FF D h / a (z) once integrated over the transverse mo- 
mentum. The Gaussian widths (k\) = 0.25 GeV 2 and (p 2 T ) = 0.20 GeV 2 [17]. The Sivers function f^(x,k\) in 
SIDIS process will be parametrized as 

}^{x,kl) = -M q (x)h(k ± )f q/A (x,k 2 x ), (16) 
where the extra fcj_-dependence h(k±) is given by 

h(kl _) = V2-e^-e-^, (17) 

with M the nucleon mass, and Mi a fitting parameter. The x-dependent part M q {x) will be parametrized as 



n (n 4- ft )( a i+i 3 i) 
M q {x) = N q x a '(l - x)^ K q+ 2[ B (1 - r, q x) 



(18) 



Compared with the previous SIDIS fits in Refs. [Tt], E^K the new ingredient lies in the factor (1 — T) q x), which is 
inspired from DSSV global fitting for the helicity PDFs [5ll. [52| • This is a simplest form that can allow a node in the 
x space: if r\ q > 1, we will have a node for x € [0, 1]; on the contrary, it rj q < 1, then no node in the region x £ [0, 1]. 
In our fit, to satisfy the positivity bound for Sivers function, we have to require < 1. To achieve this, we make 

the following substitution: 

N q ^N q /max{l,\l- Vq \}, (19) 

in Eq. (fTE)) and allow N q to vary only inside the range [—1,1], this enforces the positivity of Sivers function in x £ [0, 1]. 

Now through the relation between T qy p(x, x) and the Sivers function fy^{x, kj_) in Eq. (JTJ), we could thus obtain a 
parameterized form for T qt p(x, x) as 

T q A*,*) = ,; |2 V [jM x)f q/A {x). (20) 
({hi) + M() 

In other words, once a parametrization of the quark Sivers function is given, we automatically have a parametrized 
form for ETQS matrix element T q ^p(x, x). With this in hand, one will be able to make a simultaneous fit of both 
SIDIS at low Ph± data and pp inclusive hadron production at high Phx- As a first attempt, we will only consider u 
and d quark flavors, and include only pion data i^' ) in our fit. Thus we have N q , a q , j3 q , and r\ q for both u and d 
quarks, and M\ as our parameters, in total 9 parameters, to be determined by fitting the experimental data. 
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TABLE I. Best values of the free parameters for the u and d Sivers distribution functions, as obtained by simultaneously fitting 
HERMES [H and COMPASS data on the A^ 0h-<M asymmetry and the STAR ;53j and BRAHMS 0] data on the A N 
asymmetry. 





X 2 /d.o.f. 


= 3.6 






= 1 


N d 


= -1 


a u 


= 0.8 


ad 


= 0.8 




= 1.5 


Pd 


= 1 


Vu 


= 2.8 


Vd 


= 


Ml 


= 0.7 GeV 2 







B. Description of the data and discussion 



As we have emphasized in last subsection, the TMD factorization formula of Eq. ((6]) will be used to describe 
HERMES experiment proton target data [46J and COMPASS experiments data [47J on deuteron target. The twist-3 
factorization formula Eqs. (|9|ll|12|l will be used for An data from STAR [53J and BRAHMS [54f experiments. We use 
GRV98LO for the unpolarized PDFs [55[ , and DSS parametrization for the unpolarized FFs 56] . In our theoretical 
formalism we choose the factorization scale to be equal to the renormalization scale: n = Q for SIDIS and /j, = Ph± 
for proton-proton data. 

The results we obtain for the 9 free parameters by fitting simultaneously HERMES and COMPASS data sets on the 
Sivers asymmetry A s ^^ h and the STAR and BRAHMS data sets on the SSAs An for both charged and neutral 
pions, are presented in Table HI The extracted first moments of the Sivers functions for both u and d quark flavors 
are plotted in Fig. [TJ One can see that rjd = thus no node for d quark Sivers function. While r\ u — 2.8 > 1, there 
is a node for u-quark Sivers function that is located at Xnodc = 0.36. This value is at the border of the region probed 
by SIDIS experiments and thus in principle the node cannot be excluded. N u — 1, N c i = — 1 and thus the positivity 
bound is saturated by these functions at least at some point of x valuc0. 




-0.06 




-0.02 




(a) (b) 

FIG. 1. First moment of (a) u quark Sivers function and (b) d quark Sivers function as found using parameters of Tabled) 
here f^f m (x) = -T q>F (x,x)/2M. 



Let us now turn to the actual description of the experimental data. Our rather large x 2 with x 2 /d.o.f. = 3.6 has 
already indicated that the overall quality of our fitting is quite poor. In Fig. [2] the result of the fit compared to ir + 
HERMES [ifl and COMPASS [42} data as a function of xb- For ir~ asymmetry and Zh and P^j_ dependences the 
description is comparable to that of Fig. [2] In other words, our description of SIDIS data is satisfactory, \ 2 /#data ~ 
1.5. 



2 Since N u = 1 and 1 1 — rj u \ > 1, we have |A/"«(a:)| < 1 in Eqs. II18I19H for the whole x region except for x = 1 where the bound is saturated, 
but both Sivers function and the unpolarized PDF vanish, and thus the positivity bound is satisfied for M-quark Sivers function. On the 
other hand, JV^ = —1 and at the same time jj^ = 0, checking Afd{x) in Eqs. (I18I19D we find that the positivity bound could be saturated 
at some point of x 7^ for d-quark Sivers function. 
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(a) (b) 
FIG. 2. Description of (a) HERMES [H and (b) COMPASS [H data on tt+ production as a function of x B . 

In Fig. |31 we compare the fit with the STAR tt° data as a function of Xp for y/S = 200 GeV at y = 3.3 (a) and 
y = 3.7 (b), respectively. The solid curves correspond to the scale fj, = Phj_- The description is reasonably good, 
though slightly worse than those for SIDIS data. We have also explored the theoretical uncertainty coming from 
the scale fj, through its variation by a factor 2 up and down relative to the default values, and they are plotted as 
dashed and dotted curves in Fig. [3l This uncertainty is indeed quite large as one might expect since we are using 
the leading order formalism. The improvement could be achieved once the next-to-leading order calculations are 
performed [57j - [6ll | . 




0.1 U ' 1 1 1 1 1 1 — 1 

0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 

(a) (b) 

FIG. 3. Description of STAR n° data [HJ at rapidity (a) y = 3.3 and (b) y = 3.7 at y/S = 200 GeV. Solid curves correspond 
to the scale fj, = Ph±, while dashed and dotted ones correspond to /j, — Ph±/2 and fi = 2Ph±, respectively. 

In Fig. 21 we compare the fit with the BRAHMS tt+ and n~ data at forward angle 9 = 4° at y/S = 200 GeV. 
It is clear from these figures that our fitted parametrization for the Sivers functions (or the ETQS functions) is not 
compatible with the BRAHMS ir + and 7r~ data, and even the signs for are opposite. It is worth poin ting out that 
the previous measurements for the charged pion production in pp collisions (e.g., those from E704 [62l. l63l|) do have 
the consistent signs with BRAHMS. This finding is consistent with the heart of the "sign mismatch" paper 25]. Our 
starting point for the possibility of node in x is based on the fact that SIDIS and STAR data probe slightly different 
x region: x — xb < 0.3 for SIDIS, while x > 0.3 in the integration for Xf ^ 0.3 for STAR data. Thus a node in x can 
describe both SIDIS and STAR data rather well. However, the BRAHMS has a x region [0.15,0.3] which 

overlaps with the SIDIS data. Thus the node in x can not be a solution for the "sign mismatch" problem. Our failed 
attempt of the global fitting of both SIDIS and proton-proton data has once more confirmed that these data are not 
compatible with each other, if we consider only the twist-three contributions from the polarized nucleon. It indicates 
that there should be a sizable and more important contribution from the twist-three fragmentation in the produced 
hadrons Q- 



-0.1^ ' ' ' ' ' ' ' — 

0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 



-0.05 



-0.1 



-0.05 



-0.1 



t + + 



0.2 



0.25 



0.3 



0.35 



0.2 0.25 0.3 



0.35 



(a) 



(b) 



FIG. 4. Description of BRAHMS tt+ (a) and tt" (b) data [H| at forward angle 9 = 4° at \/S = 200 GeV. Solid curves 
correspond to the scale [i = Ph± , and dashed and dotted ones correspond to fi — Ph± /2 and fi = 2Ph± , respectively. 



Even though wc fail to cure the "sign mismatch" problem from the node in x scenario, the concept that the Sivers 
function does not need to have the same sign in the whole kinematic region (either x or k±) has important implications, 
especially when it comes to check experimentally the sign change of the Sivers function from SIDIS to DY processes. 
In Fig. [5j we show the calculation of DY asymmetry for RHIC kinematics at y/S = 200 GeV as a function of Xp. The 
solid curve corresponds to the calculations using the Sivers function with a node in x from Table. U and the dashed 
curve is the calculation based on the Sivers function from [19| which has no node in x. One can see that the prediction 
changes drastically in case node is present, however in the region of < xp < 0.25 the sign of the asymmetry is 
consistent and dictated by the Sivers function constrained from SIDIS measurements. Regardless of possible nodes 
this region is safe for measurement. In the future DY experiments, the Q 2 range will also be different, to have a solid 
prediction, one of course also needs to include the effect of the evolution [3l], [64| . 




FIG. 5. Prediction of Drell-Yan asymmetry for RHIC kinematics p'p — > £ £~ X, < y < 3. Solid line corresponds to Sivers 
function with a node from this work and dashed line to Sivers function without node from Ref. fl9t ] . The same convention for 
the hadronic frame and asymmetry is used as in Ref. [30l | . 



C. Exploration of node in k±: the simplest study 



The Sivers function with a node in k± has also been suggested as a solution to the "sign mismatch" problem in 
Ref. The main idea comes from the fact that HERMES and COMPASS SIDIS data are mostly relevant for the 
extraction of the Sivers function f^(x, kj_) at relatively modest Q 2 ~ 2.5 — 3.5 GeV 2 . Since the TMD factorization 
formalism is valid only for fciCQ, the data thus constrain the function and its sign only at a very low k± ~ Aqcd- 
However, to obtain a functional form for ETQS function T Qi p(x, x), one needs to integrate over the full range of k± 
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in Eq. ((TJ). Since currently we have assumed a Gaussian form for the fcj_-dependence which has the same sign for the 
whole fc_L region, the k± -integration will have the same sign of low k± part. However, if somehow the high k± region 
has opposite sign to the low k± part, this might alter the sign of the fcj_-moment in the integration and thus lead to 
the correct sign of T qt p(x, x). See Fig. [BJ for an illustration. 




FIG. 6. An illustration of the Sivers function with a node in k±. 



In this subsection we explore such a possibility. Talking about the fcj_-dependence, it is important to recall again 
that the relation in Eq. (fTJ) is subject to the ultraviolet (UV) subtraction and the adopted factorization scheme. To 
avoid such a problem, as a natural extension to the usual Gaussian form, we choose fc^-dependence as a difference 
between two Gaussian functions with slightly different widths. This is the simplest case which allows a node in k±, 
and we will explore whether this simple extension works in practice. The Sivers function will now be parametrized as 
follows: 



fir( x > k l) = -^q^HkiMq/A^.k]}, 
where we use the usual x-dependence (without node) for simplicity 



N q {x) = N q x°"(l - g )* lgl±gJ . 



a q *p q 



However, the fcj_-dependence h(k±) will be changed to (from Eq. (|17[l) 



h(h A 



2eM 



~M~i 



-kl/Af 2 2 

M 2 



(21) 



(22) 



(23) 



One has to choose M 2 > M 1; thus the \ow-k± will be positive, i.e., follow the same sign like the usual Sivers function 
as in Eq. ([TBI) , and the fcj_-shape follows Fig. [B] 

Using Eq. (fTJ), we could derive the functional form for the ETQS function, we have 



T q<F (x,x) = \/2e(kx) 



Ml 



Ml 



[{{kl) + Mlf 



In order that the sign of T q ^p(x,x) is altered, we should have 



Ml 



({k 



> 



Ml 



Mir 



N q {x)f q/A {x). 



We could also derive the expression for the Sivers asymmetry in SIDIS process, and find 



7 sin(tp h -4> s ) 
UT,T 



X B Z h Ph± ^ e 2 q Nq (x B )fq/A {x B )D h/q ( 



Zh 



2e(k 2 ± ) 



(24) 



(25) 
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where the widths (Pu±i) (* — 1) 2) is defined as 

(kl)M, 



(PLi) = {p 2 T)+4 { ^'^ r (27) 

In order that the asymmetry follows the same sign at low Ph± like before, one requires 

Alj Ml 
(Af| + 7 (z)(fci)) 2 < {Ml +1 (z)(kl)f' 

where 7(z)~ 1 = 1 + z,j(fc^)/(py). Thus we have three requirements from Eqs. (1251) . (|28l) plus M2 > Mi. One also 
need to take into account the fact that the Sivers asymmetries measured by both HERMES and COMPASS do not 
change sign up to Ph± ~ 1 GeV. All these requirements have constrained the allowed parameter space for Mi and 
M 2 to a very limited (small) region. For an illustration, see Fig. [7] for a typical Ph± = 0.5 GeV and z% = 0.5. This 
region gets even smaller if Ph± increases and/or Zh decreases. From such a simple study, we find that our simplest 
extension to allow a node in k± seems not to be a natural solution to the "sign puzzle" . Of course, other types of k± 
dependence which also has node in k± might still be possibl^H At the end, we want to emphasize again that there is 
the important UV regularization issue, which is out of the scope of our current study. 






0.0 0.5 1.0 



FIG. 7. The allowed region for parameters Mi and Mi at Ph± = 0.5 GeV and zh = 0.5. 

To finish this section, we study the possibility of a node in x or a node in k± in the Sivers function. The simplest 
extensions to contain such a node seem not to be working for both cases. This strongly suggests that there could be 
a sizable contribution from the twist-three fragmentation in the single inclusive hadron production 7], and we hope 
future experimental study could give us clear answers. 

IV. CONCLUSIONS 

In this paper, for the first time, we make an attempt for a global fitting of both SIDIS and proton-proton data on 
the spin asymmetries. We use a TMD factorization formalism to describe the SIDIS Sivers asymmetry for the hadron 
production at low Ph± , and the collinear twist-three factorization formalism for the proton-proton data on the single 
inclusive hadron production at high Ph± ■ We adopt a more flexible functional form for the Sivers function in order to 
simultaneously describe both SIDIS and proton-proton data. By including only the contribution from the so-called 
soft gluon pole ETQS function in the polarized nucleon, we find that we are not able to describe well all the data. 
While all SIDIS data and STAR data on tt° production can be explained, node in x does not account for BRAHMS 
proton-proton data for the 7r + and 7r~ production. We have explored the possibility of a node in x or a node in k± for 



3 We have also explored another fc^-dependent form |65|l : h(k±) = \/2e-^-e k ±/ M T (1 — r?k^) and wc find that the allowed parameter 
space for Mi and r\ is again very small. 
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the Sivers function. Our simplest extensions to allow a node for both cases seem not able to cure the "sign mismatch" 
problem. If leaving behind the UV regularization issue for the relation between TMDs and collinear functions, we 
will conclude that there could be a sizable contribution from the twist-three fragmentation function for the single 
inclusive hadron production in proton-proton collisions. We hope the future experiments could give us clear answers. 

A side-effect (lesson) learned from our study is that one should be very careful in extrapolating the Sivers function 
(or any other TMDs which do not need to be positive) to the region where there is no experimental measurements 
or constrains. The size and the sign of the functions in these region should be carefully measured in the future 
experiments. For example, careful analysis of SIDIS data at large values of xb is needed in order to rule out or 
confirm the possible node of the Sivers function. Also measurement of ir^ An at larger values of x p is needed in 
order to confirm that the node (in x) is not compatible with BRAHMS data. Careful analysis of k± dependence of 
TMDs is also needed. 
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